Rotation in an asymmetric multidimensional periodic potential 

due to colored noise 

Avik W. Ghosh* and Sanjay V. Khare''' 
Department of Physics, Ohio State University, 174 West 18th Avenue, Columbus, OH 43210 



Abstract 

We analyze the motion of an overdamped classical particle in a multidimen- 
sional periodic potential, driven by a weak external noise. We demonstrate 
that in steady-state, the presence of temporal correlations in the noise and 
spatial asymmetry within a period of the potential could lead to particle ro- 
tation. The rotation is a direct consequence of a change in sign of the noise- 
induced drift motion in each dimension. By choosing different potentials, we 
can generate a variety of flow patterns from laminar drifts to rotations. 
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Recently there has been intense activity in the analysis of stochastically driven ratchets 
These ratchets are spatially periodic systems where a spatial asymmetry in the poten- 
tial imposes a directionality, while "memory" effects from a temporally-correlated stochastic 
force (colored noise) or explicit time-dependences in the potential itself break detailed bal- 
ance, thereby rectifying microscopic fluctuations to generate a unidirectional particle drift. 
This type of unidirectional drift has been verified through experiments on colloidal particles 
or polystyrene spheres |]9|JlO[|, cold rubidium atoms in an asymmetric optical lattice [|ll| and 



quantum dots [|T^]. Predictions have also been made for SQUIDs [ITSf. The applications of 



these concepts have been manifold: this mechanism has been proposed as a possible expla- 
nation for the long-range cellular transport of motor proteins, utilizing temporal fluctuations 
at a microscopic scale to effectively transduce chemical energy into directed mechanical work 



1^ . In addition, such ideas are of interest in the nanoscale fabrication of devices. 

Most recent analyses of ratchets have been limited to one- dimensional (or "linear" ) molec- 
ular motors |T^. In this paper we show that for a particle in a multidimensional potential, 
a time-correlated noise can break detailed balance and generate rotations. Beyond their 
intrinsic interest and novelty, such rotations could be used to design non-equilibrium molec- 
ular engines. To demonstrate rotations, we develop a Fokker-Planck equation in arbitrary 
dimensions for a weak Gaussian random noise and explicitly demonstrate the existence of 
rotation in the presence of spatial asymmetry and temporal correlations. Our argument can 
be simply described as follows: in the presence of spatial asymmetry in one- dimension, a 
time-correlated noise is known to produce a drift the direction of the drift is deter- 

mined by the sense of the potential asymmetry. For a potential in multiple dimensions, the 
sense of the potential asymmetry along one coordinate can be reversed by varying the other 
coordinates, leading to a change in sign (or reversal of asymmetry) of the potential. As 
indicated in Fig. |I], the coordinate-dependent reversal of the one-dimensional drifts could 
then conspire together to generate a rotation. Contrary to rotation generated in a potential, 
which is determined by the initial conditions of the particle, the sense of our rotation is 
given entirely by the potential asymmetry, and is in fact independent of any initial condi- 



tions. Furthermore, removing either the potential asymmetry or the correlation in the noise 
destroys the rotation. 

We put our general argument outlined above in a mathematical form as follows. We 
start with a Langevin equation describing an overdamped particle driven by a Gaussian 
distributed colored noise of arbitrary correlation function in an arbitrary multi-dimensional 
periodic potential. For weak noise and small correlation time, we derive a Fokker-Planck 
equation for the probability density of the particle. We then show that in steady-state, the 
curl of the current density cannot be zero everywhere in space as long as the correlation time 
is non-zero and the periodic potential is asymmetric. Finally we demonstrate that relaxing 
either of these two conditions leads to zero current density. 

The Langevin equation describing the dynamics of an overdamped particle in a multi- 
dimensional space r in the presence of a potential U{r) and a noise f{t) is given by ||15|| : 

-fit) = w{f) + m (1) 

where dot represents a time-derivative and W = —VU (f) is the force exerted by the potential 
U{f). The noise is assumed to be Gaussian distributed, with a correlation time rf and 
strength Di along the ith coordinate. Generalizing Fox's functional calculus approach []TB 



we express the probability distribution P{f) in terms of a functional integral over different 
realizations of the noise with a Gaussian distribution. The noise in a particular direction 
i is temporally correlated through a correlation function Cj, while being uncorrelated with 
other directions j ^ i'- 



P[/] = iVexp ' ' -^-'^ /*(^)/^(^') 



— J ds J ds'y^^ 



{h) = 0; (mm) = f^-^V^r (2) 



is a normalization constant for P, and the subscripts i, j correspond to different spatial 
components of the multi-dimensional vectors. 
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The above representation of the probabihty distribution allows us to write down the 
equation of motion for the probability distribution P{r) which implicitly depends on time 
through its coordinate and the Langevin equation (|1|). The equation takes the form of a 
continuity equation P = —V • J. In the weak noise limit, the correlation time rf is much 
smaller than the diffusion time = L^/ Di over one period Lj of the potential, and the time 
t7 = LI/Uq for an overdamped particle to fall from a potential Uq. Under these conditions, 
the current density J satisfies the Fokker-Planck form [ITTD : 



f c^2 



1 + fi[T[Mii - H^/x^R - M^),, + 0(r: 



c\3 



(3) 



where the matrix elements Mjj = dWi/drj, Rjj = Y.k^kd'^^i/^rjdrk, and fi\ and /Xg are 
the first and second moments respectively of the correlation function Cj Ijl^. The role of 
color is thus to make the effective diffusion coefficient Bj position-dependent in a well-defined 
manner (i.e., a function only on the potential U{r)). 

We solve Eq. (j^) for J in steady-state (P = 0) and impose periodicity on the probability 
density P. In conjunction with the periodicity for U, this gives us a set of integral equations 
for J. In particular in two dimensions, these read: 







dxUx, y)e-'^^("'^) = P(0, y)e,(0, y) [1 - e'-^^^^-^)] 
' dyJy{x, y)e-'^y^'''y^ = P{x, 0)Qy{x, 0) [1 - e''^^^^'^^)] 
4>x{x,y) = f dz^""^^'^^ 



lo <dx{z,y) 

fy Wy{X,Z) 

(py{x,y) = / dz-f—— (4) 



Qy{x, Z) 

where L^^y are the periods along x and y directions respectively, and the coordinate zero is 
an arbitrary reference point on the x ~ y plane. 

Simplifying Eq. (^) elucidates the role of the correlation and the asymmetry terms, 
which sit on the right hand side. Expanding the square brackets to the right to the first 
significant order in rf (rf = r^, r| = r^) yields [1 — e^*^] ~ 0, where is given by: 



6,(L,,l/) = -^ / dxW,{x,y) 

Urr- Jo 



[4 (dWA' 
2 [dy 



(5) 



V 4 2 J dx^ 2 ^ dxdy 
and an analogous equation for Ly). This term is non- vanishing as long as the noise is 
correlated (rf ^ 0) and the potential in two dimensions is asymmetric, i.e., the integral over 
one period (Eq. ^) is non- zero even though the integrand itself is periodic [|T^] . 

The origin of rotation can now be traced to the existence of color and asymmetry, which 
makes the integral in (^) and thus the right-hand side of Eq. non-zero. This prevents 
J from being identically zero or constant. Now, in steady-state the divergence of J is zero, 
so the only way for J to not be identically zero within periodic boundary conditions is to 
have the curl of J not identically zero. This necessitates the current density to have local 
rotational fluxes. On the other hand, if we take the limit of white noise (rf = 0, V i) or 
make the integral in Eq. (|^) vanish (symmetric potential), the right-hand side of Eq. 
is zero, and then J can have zero curl. In fact, we can make the statement stronger: for 
white noise, the effective diffusion constant Gj = Di is position-independent. In that case, 
the Fokker-Planck equation and the steady-state lead to the following set of equations: 

V ■ J = 

V X (Je^(^^)/^) = 0. (6) 



These equations, with periodic boundary conditions and detailed balance, lead to J = [T^ 



which corresponds to an equilibrium Maxwell-Boltzmann distribution for the probability 
density P ~ exp [—U / D\. In other words, in white noise in steady-state, there is no current 
density whatsoever. The damping in Eq. (|l]) breaks time-reversal symmetry in the system, 
while color breaks detailed balance, producing rotations. 

Having thus established the necessary existence of rotational fluxes in a colored noise, 
we now move on to a concrete example for an asymmetric potential in two dimensions as 
shown in Fig. ^ The parameter a is the measure of the potential asymmetry. This potential 
is identical with respect to the transformation x ^ y. We choose our correlation function 



for both the x and y variables to be decaying exponentials in time, and assume in addition 

= t'^ and Di = D to remove any superficial differences between the x and y directions. 

We fix the boundary conditions of our potential by fixing P(0, y) = P{x, 0) = constant 
for the sake of definiteness. For a particular instantiation of our general arguments for non- 
zero rotation, we make an ansatz at this stage: we assume that is a function of y only 
and Jy is a function of x. This ansatz is consistent with the form of Eq. @ and trivially 
satisfies the steady-state condition V ■ J = 0. Then Jx{y) is proportional to (f)x{Lx,y) of 
Eq. (H) divided by J^'' dx exp[—(f)x{x,y)], with an analogous equation for Jy{x). For the 
potential in Fig. |, Eq. (|) simplifies to (p^iL^.y) = -3(r'')^f/oa7r[4sin2y sin4?/]/4i:). 

Figure |^ shows the resulting fieldplot for the current density J over one period of the 
potential. One immediately sees local rotational fiuxes separated by saddle points. There are 
drifts along x for fixed y which change sign as described schematically in Fig. [l| and produce 
a local rotational fiux. The role of spatial asymmetry and temporal correlation is clear by 
inspection of the expression for (p^iLx-, y)] the current density J is zero if we set either or 
a to be zero. In other words, the correlated noise breaks detailed balance, thereby exploiting 
the spatial asymmetries in the potential to produce local drifts and rotations. 

The relation between the potential profile in Fig. ^ and the fiow pattern in Fig. ^ can 
be summarized as follows: for a fixed x (y) coordinate, the sense of the drift along y (x) 
is given by the asymmetry integral in Eq. (^. Note that the rotations take the particle in 
and out of potential hills and valleys, implying that energy is not conserved in the rotation 
process. The particle does not follow any of the equipotentials, neither does the direction of 
the rotation depend on the initial conditions of the system. The noise-driven rotation is a 
nonequilibrium process which disobeys the fiuctuation-dissipation theorem, and the sense of 
the rotation is determined entirely by the moments of the correlation function and potential 
asymmetries. 

Varying the potential leads to a whole range of steady-state fiow patterns, including lam- 
inar fiow (global drift), global rotation (over one unit cell of the potential), and their various 
combinations. To get a drift in any one coordinate, for example, we need an asymmetric 
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potential (such as a ratchet) in that direction. Combining asymmetric potentials in various 
ways, we can generate the following classes of noise-driven fluxes [0 : 



Flow pattern Potential Form 



Rotation Ratchets coupled in x and y 

Laminar flow Decoupled ratchets in x and y 

Rotation + net drift Coupled ratchets asymmetric under 

X ^ —X, y ^ —y 

Our arguments are generalizable to higher dimensions; breaking other kinds of symmetries 
in A^-dimensions will lead to a complex network of flow patterns. 

The broken detailed balance and spatial asymmetry essential for current generation could 
be introduced in a variety of ways; in our case, we consider correlated noise and asymmetry 
in the potential. One could obtain similar results by introducing both spatial asymmetry 
and temporal correlation in the potential itself, leading to currents in a "flashing potential" 
0. On the other hand, one could include both of these in the noise itself and leave the 
potential symmetric with respect to x and y individually. This can be accomplished, for 
example, by a noise that has more kicks on average in one direction than in the other, 
thereby incorporating both broken spatial symmetry and detailed balance |]19| . 



One-dimensional ratchets have been explored under a variety of experimental conditions 
p| . To observe noise-induced rotation in higher dimensions, we now propose a realistic 



experimental set-up. Consider a system of 0.07 - 0.1 /xm charged, fluorescent polystyrene 
beads suspended in an aqueous solution at room temperature in a two-dimensional potential. 
The periods of the potential are chosen to he = Ly ^ 1 fim, with an asymmetry fraction 
~ 0.4 and maximum energy Uq ~ 75 meV. A set of crisscrossing electrodes is lithographically 
patterned to generate the two-dimensional ratchet potential similar to the one-dimensional 
potential in Alternatively we can build a ratchet optically as in Ref. [0, for example, 
by passing a low intensity laser light through a patterned reticle. For a colored noise with a 
10-40 Hz bandwidth generated electronically or optically, the particles will settle into slow 
circular orbits of ~ 1 fim diameter. Our calculations give us an estimated period of rotation 
of 1-10 hours [1^. Such fluorescent vortex patterns should be observable using a microscope. 

We have analyzed the steady-state dynamics of an overdamped classical particle in an ar- 
bitrary multi-dimensional potential driven by a noise with an arbitrary correlation function. 
For non-zero temporal correlations and asymmetries in the potential, current production 
in terms of rotations and drifts is expected. We have demonstrated rotation explicitly in 
two dimensions in the limit of small correlation times. Such rotations could be prototypes 
of periodic nonequilibrium processes such as molecular engines. By suitably tailoring the 
potential, one can generate a host of nonequilibrium flow patterns of the particle, including 
global drifts and rotations in combination. We have proposed a way of monitoring rotations 
in a realistic experimental set-up. 
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FIGURES 

FIG. 1. Schematic description of rotation over one unit cell of a two-dimensional periodic 
potential, caused by spatial asymmetry and temporal correlations. One-dimensional drifts are 
produced by asymmetric potentials in x and y, in conjunction with correlated noise. The drifts 
switch directions owing to coordinate-dependent changes in overall sign of the potential asymmetry, 
and together produce rotation. 

FIG. 2. Two-dimensional contour plot of a potential with translational asymmetry in x and y. 
The potential chosen is given by U{x,y) = Uo[smxsmy — a sin 2j; sin 2y] with Uq = a = 1. Dark 
areas correspond to valleys and bright patches correspond to hills in the potential landscape. 

FIG. 3. Local rotations (J) produced by drifts arising out of translational asymmetries in the 
potential in Fig. |2| and correlation in the noise. There is no global rotation, and the steady-state 
flow consists only of cycles and saddle points. At fixed values of x there are global drifts in y that 
reverse sign so as to generate the rotations. The arrow lengths are proportional to (r'^)^C/ga/L), so 
that as a or r"^ ^ 0, the arrow length shrinks to zero and the local current densities vanish. 
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